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Abstract 

The worldline casting of a gauge field system with spin- 1/2 matter fields has pro- 
vided a, particle-based, first quantization formalism in the framework of which the 
Bern-Kosower algorithms for efficient computations in QCD acquire a simple inter- 
pretation. This paper extends the scope of applicability of the worldline scheme so 
as to include open fermionic paths. Specific algorithms are established which address 
themselves to the fermionic propagator and which are directly applicable to any other 
process involving external fermionic states. It is also demonstrated that in this frame- 
work the sole agent of dynamics operating in the system is the Wilson line (loop) 
operator, which makes a natural entrance in the worldline action; everything else is 
associated with geometrical properties of particle propagation, of which the most im- 
portant component is Polyakov's spin factor. 



1. Introduction 



Consider a gauge theory with spin- 1/2 matter fields, in the defining representation of the 
symmetry group, described by a standard (renormalizable) Lagrangian density. A worldline 
transcription of such a system is based, _after a Wick rotation to Euclidean space-time, on a re- 
casting of the form / Vtf(x) Vip(x) e s(tf(*).tf(*M/.(*)]... — ► JVx(t) Vp(r) e s[»W.pW.^(^W)l... 
according to which the matter fields register as point-like objects whose propagation (in 
phase-space) is furnished by a path integral. To the extent that the functional integration 
over the fermionic fields is Gaussian (in the Grassmann sense), nothing is lost, in principle, 
from the content of the original system through the above transcription. Pioneering work 
that dealt systematically with the subject appears in Refs [1-3]. 

More recent advancements in worldline methodology, especially in connection to gauge 
theories, have been made, independently, by Fried [4-6], Strassler [7] and our group [8- 
12]. Respective methodological tools adopted by the above authors are: (a) Schwinger's 
functional calculus formalism [13], (b) the string- inspired computational algorithms of Bern 
and Kosower [14] and (c) Polyakov's description of the propagation of particle-like excitations 
in (Euclidean) space-times [15]. Concerning specific objectives and/or applications aspired 
to by these attempts, the following rough picture emerges. Fried and collaborators [4-6] 
as well as the present authors and collaborators [10-12] have so far focused their efforts on 
eikonal-type descriptions in gauge field theories (with QCD as the ultimate target), along 
with associated issues such as scattering amplitudes at asymptotically large collision energies, 
factorization of soft physics, heavy quark effective theories, etc. In Strassler's case, on the 
other hand, emphasis from the beginning was placed on the development of new, more 
efficient rules for performing loop calculations in pQCD as per the realizations of Ref [14]. 

A key aspect of the Polyakov-inspired approach, that we have been advocating, is its 
ability to accentuate geometrical features of particle-like propagation in space-time. To 
illustrate this occurrence consider, first, the direct performance of the path integration over 
Pp(r)R. Adhering to the use of the first order form for the Dirac operator we were able 
[9] to reproduce Polyakov's spin factor which characterizes the closed path propagation of a 
spin-1/2 particle-like entity in Euclidean space-time. This quantity has a purely geometrical 
content (it is, in some sense, associated with torsion) and accounts in a self-consistent way, 
for the re-entry of a closed (Euclidean) space-time contour by the spin-1/2 particle. By 
contrast, all other approaches to worldline descriptions have adopted the second order form 
for the Dirac operator which forces a, conventional type, spin-dependent term of the form 
<7 llv F ilv in the worldline action. Our first undertaking, in this paper, is to establish a direct 
connection between Polyakov's geometric spin-factor and the "dynamical" a ■ F term. 

A second consequence resulting from our approach is its adaptation to configurations 
involving open worldline contours. In physical terms this means that we are in a position 
to extend our considerations to processes wherein spin-1/2 entities enter as external particle 
states^. By contrast, Strassler's applications pertain to n-point functions with gauge field 

^n fact, the characterization "worldline" strictly applies to the resulting structure which, through this 
integration, formulates itself in coordinate space. 

2 It is in this context that common grounds have been established with the eikonal-based work of Fried 
[4-6], who employs Schwinger's functional calculus as methodological tool and explicitly accommodates spin 
through the a ■ F term. 
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modes as external states which attach themselves to a fermionic loop [7]. The major part 
of our present effort is to extend worldline computational algorithms, hitherto provided 
for diagrams involving fermionic loops, to ones wherein open fermionic paths make their 
entrance. 

Our ultimate hope is that, through the present work, the elegant separation between 
purely geometrical properties of matter particle paths, on the one hand, and dynamics, on 
the other, will emerge as one of the most potent aspects of the worldline casting of gauge field 
theories. Leaving the elucidation of this remark to an a posteriori discussion, let us proceed 
to exhibit the organization of the present paper. In the next section we display worldline 
expressions which will form the basis of our subsequent considerations. We shall proceed, 
in Section 3, to make the connection between the "dynamical" a • F term, appearing in the 
worldline action as a result of adopting the second order form for the Dirac operator and 
the "geometrical" Polyakov's spin factor. This is accomplished via a fortuitous application 
of an area derivative operator acting on functionals [16,17]. Section 4 constitutes a central 
part of our paper. Specifically, we undertake the task of carrying out the procedure that will 
lead to an explicit worldline formula for the fermionic Green's function to second order in its 
perturbative expansion. A major aspect of this attempt is to present an expression which 
employs super-particle coordinates, i.e. a path integral description which uses both bosonic 
and fermionic variables. We proceed, in Section 5, to discuss computational tools which are 
involved in practical applications of the worldlike scheme. A first such undertaking is carried 
out in Section 6 where we establish, to second order, the equivalence between the Feynman 
diagrammatic and worldline perturbative expansions of the fermionic propagator in a non- 
abelian gauge theory. A direct one loop calculation is performed in Section 7; it produces a 
final expression in parametrized form which is ready to enter perturbative estimates of higher 
order Green's functions with external fermion (as well as gauge field) legs. Our concluding 
remarks are presented in Section 8, while some technical matters pertaining to the work in 
Section 7 are displayed in an Appendix. 

2. Basic worldline expressions 

In this section we shall discuss worldline formulas placing emphasis on fermionic Green's 
functions in a background of, non-abelian in general, gauge fields. We shall work in Euclidean 
space-time in the context of which the spin-factor can be defined. Special emphasis will be 
placed on the second order Dirac formalism as our immediate aim (cf. next section) is to 
establish its connection with the spin-factor. 

Let us commence our considerations with a short presentation of previous results stem- 
ming from the employment of the first order Dirac operator. A systematic procedure which 
takes one from the original field theoretical casting of the system to worldline formulas can be 
found in [9]. Here we shall proceed formally. Consider, first, the effective action functional, 
with dependence on the background gauge field. Its formal expression is 

W[A] = lnDet[i-y ■ D - (m - ie)], D„ = d„ + igA^. (1) 
Using Schwinger's proper time representation [18] and going to Euclidean space-time we 
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write 

W\A] = -Tr I" tEe-n-TD+m) (2) 
Jo T 

Computation of the non-Dirac 7 matrix as well as non-color part of the trace in coordinate 
space representation leads, after a suitable choice of gaugeQ, to the worldline result [9] 



poo fi'i • /• pT 

W\A] = - / — e~ Tm / Vx(t)6[x(t) 2 - l]tr$(C)Tr c Pexp ig x- A(x(t)) 

JO T Jx{0)=x(T) |_ JO 

(3) 

where tr(Tr c ) denotes trace over 7-matrix (color), P stands for path ordering and where the 
spin-factor $(C) [15] is given by 



$(C) = Pexpl- [ dr 
\ 8 jo 



unA x { T )\b^iv\\ ( 4 ) 



with 



^VOO)) = •A%p(t)x v (t) - x^{t)x v {t\ x 2 = 1, (5) 

describing the orientation tensor of the local, perpendicular to the path, plane. 

As already mentioned, our main preoccupation in this paper will be with fermionic 
Green's functions (propagators) whose definition involves open paths. To begin, let us dis- 
cuss worldline expressions for the fermionic Green's function, in a (non-abelian) gauge field 
background. We first consider the case where the first order form for the Dirac operator is 
adopted. The corresponding Schwinger proper time formula for the Green's operator is 

roo 

iG\A} = - dT e~ T{ ~ TD+m) . (6) 
Jo 

Its worldline representation, in phase-space, reads 

iG(x,y\A) = f Ve(r) f Vh{r) L m=x Vx{r) f Vp{r)exp[i [* drh(r)e(r)] 

J J J x(T)=y J JO 

xexp[i / drp(r) ■ x(t)]P exp[— / dTe{r){m + 27 • p(r)}] 



xPexp[ig A(x(r)) ■ x(r)). (7) 
Jo 

The integral over T>h(r) on the right hand side supplies our commitment to a choice of 
gauge with respect to the reparametrization invariance requirement. The latter, corresponds 
to the transformation r — » t = f(r), /(r) > 0, /(0) = 0, /(l) = 1, under which we have 
e(r) - f(r)e(t), x(r) - f(r)x(t) and h(r) - j^h(t). 

The second order formalism for the Dirac operator, on the other hand, bases itself on the 
operator expression 

iG[A] = (7 ■ D + m) I" dT e -n™ 2 -(rDf] ^ 
Jo 

3 We are here referring to a gauge choice with respect to parametrization invariance and mass identification. 
More specific comments on this issue will be made shortly. 
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whose phase-space worldline representation turns out to be 



iG(x,y\A) = I Ve(r) f Vh(r) L )=x Vx(r) f Vp(r) 

J J J x(l)=y J 



m — 



2e(l) 



7 • 



xexp 
xexp 



dr/i(r)e(r) 



xP exp 



exp[i / drpir) ■ x(t)] 
Jo 

- I rfre(r){p 2 (r) + m 2 } 
Jo 

rl If 1 

ig / A x t • x(t) + -g / dre r )a ■ F{x{t)) 
Jo I Jo 



(9) 



As already stated our first preoccupation in this paper is to establish that the a ■ F term 
in the above formula translates into the analogue of the spin-factor for open fermionic paths. 
The first step in this direction is to carry out the integration over Dp(r). 

We obtain 

1 C(e) 
j j 

1 



iG(x,y\A) 



J Ve{r) j Vh(r)exp i J dr/i(r)e(r) 

x L 0)=x Vx(T)[m - t-4tt7 • 
J x(i)= y ze(l) 



xexp 



xP exp 



1 r 1 



4 Jo 



dr 



x 2 (t) 



e r 



— m 



dreij) 

1 rl 



ig I A{x{t)) ■ x{t) + -g f dre(r)a ■ F(x(r)) 
Jo I Jo 



where 



C(e) 



n 



N 
n - 



/ N \°/ 2 
1 \47re r 



(10) 



;n) 



can be looked upon as a normalization factor. 

Our immediate objective is to recast the above expression into a form where the open-line 
analogue of the spin-factor replaces the a ■ F term. The relevant work will be carried out in 
the following section. 



3. Spin- factor in the second order formalism 



We commence our considerations which will lead to the recasting of the a -F term in (13) 
by utilizing the area derivative operator for path dependent functionals. We shall follow, 
in this respect, Polyakov's proposal [16] for practical reasons. A geometrical definition has 
been given by Migdal [17]. Let us, then, introduce the operator (area derivative) by 

X f+e S 2 

*}_ "" s (t ..,^.u (12) 



lim F 



/2)5x v (t - a/2) 



Consider the action of the above operator on the functional 



7[i;(r)] = exp 



dr- 



x 2 (t) 



(13) 
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which leads to 



exp 



4 Jo e(r) 



1 ['^(r) 



4 Jo e(r) 



with uJn U now given by 



lo^ u (x(t)) 



2e 2 (i 



[x^(t)x v (t) - Xh(t)x v (t)]. 



On the other hand, we find 



(14) 



(15) 



P exp 



1 l r\ 

ig I A(x(t)) ■ x(t) + -g dre(r)a ■ F(x(t)) 



Pexp 



/ dre(r) a 
Jo 



2 U Jo 
5 



P exp 



2 Jo Ss(t 
Inserting into (10) we readily derive, after partial integration, 

C(e) 



ig / drx{r) ■ A(x(r)) 
Jo 



i I drhMeir) 
Jo 



iG(x,y\A) = J Ve(r) J Vh{r)exp 

x L a)=x Vx{T)[m - ^jrzl ■ x(l)] 
J x (i)= y 2e{i) 

.If 1 , X 2 (t) 2 f 1 

xexpl — / dr — — m / arelr 

L 4 Jo e(r) Jo 



xP exp 



A Jo 



dr o ■ uj(x(t)) 



P exp 



ig / drx{r) ■ A(x(t)) 



(16) 



(17) 



Even though the term 'spin-factor' per se pertains, as a geometrical quantity, to closed paths, 
we shall employ, from hereon, this nomenclature for the expression given by Eq. (4) even 
when it enters quantities defined on open paths, such as the above. 

The same considerations applied to the effective action functional whose formal expres- 
sion, in the second order formalism, is 



- W[A] = l -Tr jf 
leads to the worldline representation 



dT 
o 1 



-T(m 2 -( 7 -D) 2 ) 



-W[A\ 



Ve(r) / Vh(r)exp 



i / dr/i(r)e(r) 
Jo 

l r 1 . x 2 ' 



(7(e) 



/ Vx(T)exp\ — / dr - , , m 2 / dreir) 

Jx(o)=x(i) 1 4 Jo e(r) Jo 



c(0)=x(l) 

xtrP exp 



'A Jo dTa - UJ ^ T )) 



Tr c P exp 



ig / dTx(r) ■ A(x(t)) 
Jo 



, (19) 



4 Note that the expression that follows coincides with eq. (5) when e 2 (r) = 1, which corresponds to a 
specific choice of gauge. 
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where (7(e) = ^yC(e). 

For a proper choice of gauge the above result matches the one that has been obtained 
[9] through the first order formalism, cf. Eq. (3). The notable accomplishment resulting 
from the manipulations that have just been carried out is an apparent dissociation of a 
geometrical feature (cf. spin-factor), characterizing the propagation of particle-like modes, 
from the dynamics operating in the system. The latter seem to be exclusively carried out 
by the Wilson factor P exp ig Jq drx(r) ■ A(x(t)) . As will actually turn out, computations 
involving the spin factor term do, in an implicit manner, associate themselves with dynamical 
input. 

Before proceeding further we shall commit ourselves to a choice of gauge in relation to 
reparametrization invariance. Designating as dynamical variable the proper time we set 



/ r rfr'e(r'), T= fdreir) 
Jo Jo 



(20) 



and obtain 

iG(x,y\A) 



f°° 2 f i if 

/ dTe~ Tm / B(0)=B Vx(r)[m - -T x(T)]exp -- drx 2 (r) 

JO J x(T)=y 2 L 4 JO 

xP exp — - / dTa-w(x(r)) P exp ig / drx(r) ■ A(x(t)) 
4 Jo jo 



, (21) 



where the measure Vx(r) is now defined as UmN^^ i^f) 2 Y,n=od D Xn, while w^ u is, on 
account of our choice of gauge, given by 



Note that Iwl 2 



2 , which means that it furnishes the curvature when \x\ =const. 



(22) 



Turning our attention to the full fermionic Green's function we write 

1 , 1 rl 

Vx{t) m — 

>x{T)=y 



roo 2 f I 1 

iG(x,y) = / dTe~ Tm x(0)=x Vx(t) m - -7 • x(T)]exp[-- / dTX 2 (r) 
Jo Jx(T)=v - 2 4 Jo 

r fT 

< P exp ig I drx(r) ■ A(x(t)) 
Jo 



xP exp 



dra ■ w(x(t)) 



>a,(23) 



where the expectation value in the gauge field sector includes gauge fixing terms, ghost 
integration and Dirac determinant contribution. 

Generalizing, we infer that, in the worldline language, quantities of theoretical interest, 
such as n-point Green's functions and generating functionals, draw all dynamical input from 
expectation values of Wilson lines or loops. Such objects have been extensively studied in 
the past, see, for example, [16, 19-23]. Any possible non-perturbative application of the 
present formalism will, therefore, be greatly facilitated by relevant, accumulated knowledge. 

Even though Eq. (23) will serve as the starting point for our subsequent considerations, 
we find it useful to recast it in a form which adopts particle-based coordinates. Along with 
practical consequences, this will facilitate a direct comparison with Strassler's (closed path) 
worldline expressions. We shall devote the remainder of the present section to this task. 
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As a first step, let us cast (23) in the form iG(x,y) = a(u 2 ,m 2 )u ■ 7 + b(u 2 ,m 2 ), 



±(y - x)^ by writing 



iG(x,y) 



1 r c 
D Jo 



-Tm x 



L(o)=x T>x(t) 



X 



dTe 

'x(l)=y 
m x 1 • (rr\ 



tr 



-falpP exp 



r 



dr a ■ w(x(t)) 



x < P exp 



ig / dTx(r) ■ A(x(t)) 



>a 



(24) 



which follows via the identifications b = ^tr(iG), au^ = j^tr^^iG) and once taking into 
account that trl = D. 

We now use the identity 



tr \ leTipP exp 



4 Jo 



r 



xexp 



dra ■ w(x(t)) 
1 f T 



[ 



</>(0)=V(T)=0 

- / drip ■ ip + - / rfr^(r)^(r)w^(x(r)) 
2 Jo I Jo 



[#]2^(T)^(T) 



(25) 



where the ip^ are Grasmann variables whose correlator, with respect to the action — | J T drip ■ ip 
is determined by taking into consideration that 

and 



D h 



D ^(o)+v>(T)=o 



[dip] exp 



1 f T 
-- / drip-ip 
l Jo 



We obtain 



< Vv( t i)Vv(t~2) >i,= -S^signin - r 2 ). 



(26) 
(27) 
(28) 



One is thereby led to the following particle-based, to be referred to as 'super-particle', 
representation for the full fermionic propagator 



iG(x,y) = — dTe~ m L {0)=x Vx{t) — S a/3 - -j a x (T) / [dip] 

L> Jo J x(l)=y It) 2 J Ji>(0)+i>(T)=0 



x2iP a (T)iPp(T)exp 



1 r T 



-I 

4 Jo 



drx' L 



1 

2 

1 r T 



2 Jo 



1 r T 

drip -*P + - drip^ipvW^xir)] 



2 Jo 



x < P exp 



ig I dTx(r) ■ A(xM) 
Jo 



>A ■ 



(29) 



This is an essentially new result which not only extends Strassler's expressions to open 
fermionic lines but encodes spin effects through a geometrical quantity entering a particle- 
based action functional. At the same time it underlines the role of the Wilson line operator 
as the sole agent of the dynamics operating in the system. The spin factor, in other words, 
has joined the rest of the path-dependent (super) coordinates to account for particle-based 
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characteristics of spin- 1/2 matter field propagation. Generalizations of spin factor expres- 
sions pertaining to higher spins have been discussed, from a different viewpoint, in Ref [24]. 
The bottom line is that the spin factor accounts for geometrical features of paths induced 
on them by the nature of the particle-like entity which traverses them. 

Mention should be made, at this point, of the attempt in Ref [25] to tackle the problem 
of open line propagation of spin- 1/2 particles in a gauge field background. The relevant 
methodology has a different philosophy (analytical rather than geometrical) and restricts 
itself to situation wherein the gauge field(s) do not acquire a dynamical character, i.e. they 
remain external. By staying faithful to Polyakov's geometric point of view [15], on the other 
hand, we have both attained simplicity in form of worldline expressions while achieving, at 
the same time, extensions beyond the free particle case. 

The corresponding expression for the effective action functional is 



W[A] 



2 J^(0)+V(T)=0 

r i f T 

xexp — / 
F [ 4 Jo 

x Tr c P exp 



j 00 dTe~ Tm2 L 0)=x Vx(r 

JO J^lT\=ii 



x(T)=y 
1 f T 1 f T 

drx 2 - - drip ■ ip + - drij) ^ v w ^ v {x{t)) 

_ J Z J 



ig / drx(r) ■ A(x(r)) 
Jo 



(30) 



The above equation corresponds to Strassler's result, modulo the presence of the spin- 
factor term in place of a ■ F. It serves as the starting point for computational rules involving 
quark loops. The relevant particle-based action can be read off the above formulas and has 
as follows 

S = - drx 2 + - / drip ■ 4> — - / dT%l)nip v Wn V (x{r)) — ig / drzx-Az, (31) 
4 Jo 2 Jo 2 Jo Jo 



6x^ the first 



where the z-fields serve to enforce path ordering when background fields are present 

Note that under the supersymmetric transformation 5x^ = —i8ip^, Si/j^ 
two terms in (31) remain invariant whilst the last two do not. On the other hand, we can 
easily retrieve the a ■ F term, by partial integration, from the spin-factor. This brings us 
right back to the Strassler form for the worldline action which is supersymmetric invariant^]. 
We shall refer to (31) as super-worldline particle action. 



4. Analytic manipulations with the spin factor: Fermionic prop- 
agator to 2nd order 

Besides its interpretational appeal as a geometrical agent that accounts for the spin of 
a propagating particle, the spin factor presents practical advantages as it actually expedites 
analytic manipulations involved in the computation of fermionic Green's functions and ef- 
fective action functionals. Samples of the latter case have been given in [9]. Here we focus 
our efforts on calculations pertaining to fermionic propagators. 

5 Note that supersymmetry appears to be manifestly broken in (31) due to the boundary conditions 
involved in the partial integration which produced the spin factor. 
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4a. General Considerations 



According to the resulting expressions in the previous section the spin factor expression 
combines bosonic with fermionic particle-based coordinates. Suppose we were to expand 
the relevant exponential. We would then find ourselves having to compute expectation 
values in the bosonic sector of the type < (x^{t)x v (t) — x fl (r)x u (r))... > x . Clearly, the 
only way to get a non-null result is for a four- vector, say k^, to make its entrance so that, 
in combination with the four-vector (x — y)^ which will emerge from the path integration, 
would facilitate the creation of antisymmetric combinations. Equivalently, unless the first 
derivative is discontinuous at some r = n the aforementioned expectation value terms, once 
recast the form < (x fl (r)x u (r + a) — x Ai (r)i; I/ (r — a))... > x (a — >• 0), would lead to a vanishing 
result except at points where a 'force' acts, thereby 'injecting' a four-momentum k^ on the 
fermionic line. 

Consider the situation where M such points are present on a given fermionic line. Let 

us set 

ii 

We surmise (here C denotes an open curve) 



^[x^Ti + e) - x^Ti - e)} = lim e ^ iki^ i = l,...,M. (32) 



$(C) — > lim t ^ Q exp 



M 



i=l 
M 



n+e 



T 

-J2 I dT^^T)^ v {T)[x^{T)X v {T) - XviTjX^T) 
4 JTi-e 



T 

= exp-—22i>iJ,(n)i>v(Ti)[iki tJi x l ,(T i ) - ik iu x^{Ti)}. (33) 
1 i=i 

Momentum conservation can be taken into account by setting Y^fLi hfj, — 0. 

Clearly, M measures the number of points where a gauge field is applied on a given 
fermionic pathF]. This means that M coincides with the perturbative order to which one 
expands the Wilson line in a given worldline expression such as, e.g., (29). 

Let us now return to the expansion of the exponential that furnishes the spin factor, only 
this time we focus on its fermionic component. From eq (26) we determine 

< ^(t)Vv(t) >^= -V (34) 

This means that any term in the expansion of order greater than M will inevitably lead 
to mixtures of symmetric and antisymmetric factors, hence a vanishing contribution to the 
propagator. So, not only is M associated with perturbative order but also sets a bound on 
the contributing terms in the expansion of the spin factor. 

Armed with the above concrete as well as of practical value realizations, we are in position 
to turn our attention to the 2nd order expression of the fermionic Green's function and gain 
specific insights with respect to spin factor contributions to its worldline form. This is 
precisely the task that we shall carry out in the next subsection. 



4b. Spin factor contribution to the second order fermionic Green's function 

6 It is in this context that an earlier comment, to the effect that the computation of the spin-factor 
expression requires dynamical input, was made. 
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Expanding the Wilson exponential in (29) we bring down, to second order in g 2 , the 
factor Jq dri Jq dr2X^(T 2 )x u (Ti) < A^(x{t 2 ))A v {ti)) >a- Working in the Feynman gauge we 
write 



xx(t 2 ) -i(ri) 



m K 1 • (rr\ 



[#(r)]2^(T)^(T) 

V>(0)+V(T)=0 



xexp 



\ r T 1 r T ■ 1 r T 

- I drx 2 - - drip ■ tp + - / drip^w^ - ik ■ (x(t 2 ) - x(ti)) 
4 Jo 2 Jo 2 Jo 



.(35) 



Focusing our attention on the spin-factor expression, let us consider the first order term, 
to be denoted by in its expansion. According to the general discussion of the previous 
subsection we determine 

T 

= ~z XI Vv^OVv^fe^^) (36) 

1 i=l,2 

Setting fc 2 = —k\ = k we obtain 

T 

= — Vv( r 2)Vv(T 2 )[^M r i) -*Mv( r i)] 

T 

-^-^(ti)^(ti)[^^(t 2 ) - ifc I/ i /i (r 2 )]. (37) 

Consider, next, the path integration over the ip fields. According to (34) the contribution 
from the above terms is zero. Had our computation referred to the effective action functional 
W[A] this would be the end of the story, as far as path integration over the fermionic 
coordinates is concerned. For the Green's function, on the other hand, a four-spinor factor 
also appears for which we determine (we keep only the antisymmetric contribution, under 
the exchange // <-> i/) 

2 < ^(T^^T^^T^^i) > i ,= ~[S ait Sp v - SavSfo], i = l,2. (38) 

Reserving the same notation for the quantity that will result after the path integration with 
respect to the particle coordinates we write 

T T 

-> — ik a [xp(T 2 ) -xp(n)] - —ikp[x a (r 2 ) -x a (ri)], (39) 

where the arrow serves to signify the fact that a functional integration with respect to the 
X/j,(t) fields remains to be carried out. Now, along the way, the full computation of will 
encounter the quantity x(t 2 ) • x{ji). Taking into account that 

2 

< x p (r)x a (r) > x = -5 pa , (40) 

whose justification is, as we shall see later, intimately connected with reparametrization 
invariance, we are in position to replace (39) by 

-> ik^x^n) - x p {t 2 )\ - ikf}[x a (n) - x a (r 2 )}. (41) 
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Next, we consider the contribution from the second order term which, according to our 
general discussion, is given by 

T 2 



$(2) 



8 



-^(^^(^^(n^A^i)^^ <-> 2). (42) 



Taking on the task of calculating expectation values in the ^-sector of the (super)particle 
system we determine, once keeping only the symmetric term under the exchange 1 <-> 2, 

2 < V ; a(^)V ; /3(^)^( r 2)^( r 2)V ; K(n)V ; A(Ti) ></,= ^8 aj3 [8^vK - S v \Shk] [signal - r 2 )] 2 (43) 

As before, we keep the notation $ < - 2 - ) for the quantity which results from path integration 
in the super-particle sector. Substitution of (46) into (45) furnishes the intermediate result 



$(2) 



rp2 

o 

= — -rSapi^S^ - k^k u )x^(T 1 )x u (r 2 )G 2 F (Ti,T 2 ), 



(44) 



where we have, following Strassler, set Gf(ti,t 2 ) = signiji — r 2 ). 

Further utilization of the bosonic sector path integration allows us to cast $( 2 ) in the 
more simplified form 



$ (2) - -M^<W - k^S^GUn,^) = -5 a p{D - l)k 2 G 2 F (n,r 2 ) 



(45) 



All higher order terms in the exponential expansion of the spin-factor term give vanishing 
contribution to the Green's function. Here, we explicitly witness the harmonization of the 
spin-factor with dynamical aspects of the calculation: Recognizing the order of the pertur- 
bative expansion the spin-factor series terminates when the number of vertices (points on 
the path where a 'force' is exerted) is saturated. 

Putting together everything that was so far carried out, we arrive at the following world- 
line expression for the fermionic propagator, to second order in perturbation theory, 

iG {2 \x,y) = -g 2 c F f 1 f° dTe~ Tm * f dr x f dr 2 9(r 2 - n ) L )=x Vx(r) 

J (Zir)^ Ar Jo J o Jo J x (T)= y 



x 



Tfl ! 



{S a/ 3x(r 2 ) ■ x(ti) + ikp[x a (Ti) - x a (r 2 )} - ik a [x p {Ti) - x p {t 2 )} 



-5 a(3 G 2 F (n,T 2 )(D-l)k 2 }exp 



4 Jo 



dr x 2 — ik • (x(ti) — x{r 2 )) 



(46) 



The corresponding expression for the effective action functional W[A] follows in an anal- 
ogous way and reads 

n r rjDu poo , pT pT 

W[A] = -g 2 j Tr c A,(k)A u (-k) jf dTe~ Tm2 jf d n J q dr 2 9(r 2 - n ) 



x 



lx(0)=x(T) 

xexp 



Vx(t)[x^(ti)x v (t 2 ) - G F (Ti,r 2 )(k 2 5^ - k^k v ) 
1 r T 



(47) 
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The above results express the embodiment of spin factor contribution to worldline for- 
mulas. Clearly, the immediate question is whether (46) furnishes the same perturbative 
expression for the propagator that one extracts, conventionally, through Feynman diagram- 
matic, rather than worldline, reasoning. That this is, indeed, the case will be established in 
Section 6. More interesting, of course, is the question regarding whether the formalism that 
has been developed in this work can define autonomous approaches to field theoretical com- 
putations or, at least, offer practical advantages for higher order perturbative calculations 
for processes involving external fermionic fields. Setting this issue aside we shall, in the next 
section, turn our attention to the establishment of a set of rules pertaining to expectation 
values of bosonic 'field' products. 



5. Computational tools in the super-particle representation 

The novel feature of the worldline formalism is the super-particle mode of description 
which pervades its expressions for each given quantity of physical interest. Our spin factor 
manipulations in the previous section has mainly concentrated on aspects surrounding path 
integration over fermionic coordinates. In the present section we shall conduct a systematic 
discussion of the particle-based (super) sector whose ultimate aim is to promote the practical 
side of the worldline scheme. 

Focusing on the bosonic sector of super-particle actions we realize that its kinetic term has 
the conventional form J T dTx 2 (r). Corresponding propagators can be appropriately defined. 
For open lines the Green's function is denoted by A(r, r'). It obeys the equation 

fP A(t,t')) = -5(t-t>), (48) 



dr 2 



with boundary conditions A(0, r') = A(T, r') = 0. 

For closed contours the situation demands closer inspection. As argued by Strassler [7], 
the relevant Green's function, to be denoted by Gs(r, r'), obeys an equation of the form 

^ G b(t,t')) = 5(t-t')-± (49) 

where T is the length, according to a chosen parametrization, of the (closed) contour. Fol- 
lowing Strassler we adopt the convention which sets C = 2. 
One obtains, respectively, 

A(r, /) = T -^ ¥ ^9(t> -t) + T -^ F ^9(r - r'), |U(t, r') = 9(r> - r) - T - (50) 



and 



Gb(t, t') = L^(T - |r - r'l), ^-G B (r, r>) = sign(r - r') - (51) 



We also mention the relations 
J^[A(t,t 2 ) - A(r,n)] = ~[G fl (r,n) -G b (t,t 2 )}, ^A(ry) = ^ B (r,r'), (52) 
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where the double dot on G B signifies second derivative with respect to r. Finally, during 
our calculations we shall, following Strassler, use the rules Gb(t,t) = and G B (r,r) = 0. 
G b (t,t>) = -G b (t>,t). 

As a first application let us establish the validity of Eq. (40). We readily determine 

1 f T ~ e 

< xJr)x a (r) > x = lim e ^ —5 pa dr'G B (Tr') 

1 Jt+e 

2 r T ~ e r l i 2 

= lim e ^ -S^ J dr' [£(r - r) - -J = -S pa , (53) 

One, thereby, explicitly verifies that integration over the bosonic coordinates x(t) produces a 
result induced by parametrization invariance as it yields a constant measure for the velocity: 
<|if>=f. 

Next, we undertake the computation of expectation values of the first few products of 
derivatives of bosonic 'fields' i M (r), with action functional the one that arose in connection 
with our second order fermionic propagator expression, cf. eq. (46). The particle-based 
(bosonic) action reads 

1 r T 

S = - drx 2 + ik ■ \x(t 2 ) - xin)]. (54) 
4 jo 

The fact that terms of higher than 2nd power are absent in S encourages the employment 
of techniques which base themselves on solutions of the classical equations of motion. The 
latter read 

x*(t) = 2ik,[5{r - r 2 ) - S(t - n)]. (55) 

with solution 

tf(r) = V ,(r) ~ 2i^[A(r, r 2 ) - A(r, n)], (56) 

where ^(r) = t + x p = u p t + x M . 

Under the substitution x(t) — > x(t) + x d (r) the action becomes 



1 f T 

S = Vci + 

with 



S = S cl + ] f T drx 2 (r), (57) 
4 jo 



S d = -Ju 2 + ik ■ u(t 2 - n) + fc 2 G B (n, r 2 ), (58) 
Consider, first, the case of < x^t) >. We readily determine 

< x„{t) >= AT 1 L 0)=x Vx(T)x p (T)e- s = x*(t), (59) 

J x(T)=y 

where we have defined 

M = Lo)=« ©a;(r)e- s . (60) 

J x(T)=y 



Note that 

= ?/. - Oik , 

l dr 



w M - iA; M [Gs(T,Ti) - G b (t,t 2 )\. (61) 
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which, given the equal time rules mentioned after eq (52) along with the fact that G\b(t, t') = 
-G b (t',t), yields 

*t in) = u M + ik^G B (r u r 2 ), i = 1, 2. (62) 

Observe, in passing, the consistency of the above result with substitution (32) associated with 
non-vanishing contributions of the spin factor in the bosonic sector. Indeed, we determine 

< x IJ ,{Ti + e)-x fl (T 1 -e) > x = -ik^GBiri + e!,^) -G B (r 1 -e 1 ,ri)} = —2ik fl = 2ik lfl . (63) 

For the correlator we find 

< x Ai (r 1 )x M (r 2 ) > 



where Af a = {47T ^ )D/2 ■ 

The last relation brings to surface the fact that G B (ri,Tj) will make regular appear- 
ances into our expressions. An immediate problem can be identified with its delta function 
component whose presence becomes the source of potential ultraviolet problems. The way 
to confront this issue is by performing an integration by parts which replaces GB{Ti,Tj) 
by G|( Ti,Tj). Even at that one could still anticipate problems associated with the end 
point behavior of integrations with respect to the dr^s, where Tj denotes a point of gauge 
field contact with the super-particle contour. As with the case of closed contours [7], what 
saves the situation is the fact that partial integration invariably produces the combination 
G 2 B (Ti,Tj) — G F (ri,Tj) which eliminates any pathological behaviorf]. The super-particle de- 
scription, which, in our treatment, has been introduced through the spin factor, plays a 
crucial role in producing singularity free expressions. 

In conclusion, our work in this section underlines the parallelism between open and closed 
line computational subtlelties. Nothing that arises in the former fails to do so in the latter 
case, or vice versa. The only difference is that additional terms in the super-particle action 
are induced by the new boundary conditions but which have no fundamental impact, as 
far as divergent behavior is concerned. Even though we shall not, in the present paper, 
study physical quantities in which external gauge field lines attach themselves onto matter 
particle contours we fully expect to encounter analogous computational features, e.g. pinch 
singularities, in these situations as well. We plan to report on such calculations in the near 
future. 

6. First application: Recovery of Feynman diagrammatic expan- 
sion to 2nd order 

We have already pointed out that nothing is lost to the original content of the field 
theoretical system on the account of manipulations that have taken place for the purpose of 
casting it into its worldline form. As a first illustration -which, at the same time, will enable 

7 As an explicit illustration of this occurrence the reader is referred to the next section where a computation 
pertaining to G^ 2 '(x,y) will be carried out. 



x{0)=x{T)=0 



1 N , ^2 

Vx{r)x iJL {T 1 )x l ,{T2)e'^ Jo 



tf(n)x d Jr 2 ) 



5/ivG 



b(ti,t 2 



(64) 
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us to see the corresponding particle-based computational tools at work- we intend, in this 
section, to prove the equivalence between the worldline and Feynman perturbative expansions 
for the fermionic propagator, to second order. Given the notable calculational advantages 
that have already been established for non-abelian gauge field theories in connection with 
closed particle contours [14, 7, 26], one is encouraged to think that our present extension 
will lead to efficient ways for perturbative estimates, to a given order, for processes involving 
open matter particle contours as well. 

Let us commence our considerations by dispensing with the Oth order term of the 
fermionic propagator. As per our discussion in the previous section, we expect no contribu- 
tion from the spin factor in the super particle action. Indeed, the only available four-vector 
in the free propagation case is = ^(x — y)^ and there is no way to generate antisymmetric 
combinations from it. The relevant worldline expression reads 



iG<®{x,y) = (°° dTe- Tm2 L 0)=K Vx{r)[m - U ■ x{T)]e^ C d ^. 

JO J x(T)=y * 



(65) 



The solution of the classical equations of motion is, simply, x cl {r) = rj^, with r]^ as in eq. 
(56). Under the substitution x(t) — > x cl {r) + x(t) path integration is performed with the 
same action but with boundary conditions x(0) = x(T) = 0. Given that the classical value 
of the action is \Tu 2 and in view of (59) we finally obtain 



iG^\x,y) 



e -Tm 2 y 

dT- — r—w—\ \ n /o \ m 7 



tie 



jD, 



-iq(x-y)-Tq 2 



(27T) 



D 



(m — i~f ■ q). 



(66) 



Integration over dT finally produces the free fermionic propagator. 

Focusing our attention to the mass containing term I m in iG^ 2 \x,y) we first recall that 
spin factor contributions are already contained in (46). Moreover, we have carried out, in the 
previous section, the computation of relevant expectation values with respect to the bosonic 
action given by eq. (54). Putting everything, that has already been worked out, together we 
obtain 

2 r d D k i yoo dT „ Tm2 r T , r T , Q , v 
Im = ' m9CF Jj2^¥Jo WrW e Jo dTl Jo dr ^- T ^ 



X 



{DG b (t u r 2 ) + [u + ikG B (n, r 2 )] 2 — (D — l)k 2 G 2 F (r u r 2 )} 



exp 



-ik ■ u(t 2 - n) - k 2 G B {ri, r 2 ) - 



T 



-u 



(67) 



Our immediate task is to perform an integration by parts in order to get rid of the 
undesirable quantity Gb- We readily determine 



dr x \ dr 2 9(T 2 - Tx)G b {t u r 2 )exp 



—ik 



u{r 2 - n 

T 



- fc 2 GB(ri,r 2 ) 

-T 



r ■ r i f f 

/ dr 2 G B (r 2 ,0)exp -ik ■ ut 2 - k 2 G B (T 2 ,0) + / dn dT 2 9(r 2 - n) 
Jo 1 J Jo Jo 

x[k 2 G 2 B (T 1 ,T 2 ) - ik ■ uG B (r 1 ,T 2 )]exp -ik ■ u(t 2 - t ± ) - k 2 G B (T 1 ,T 2 ) 



(6J 
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For the reader who is acquainted with closed-contour worldline manipulations the analogy 
with the above practice is familiar; simply the exponential has more terms for open lines. 
The important point is that in both situations the integration by parts goes through on the 
basis of the adopted regularization, namely Gb(t, t) = 0. 

Making the variable change r 2 — T\ = aT and r 2 + T\ = 2/3T we readily determine 



2 r d D k 1 f°° dT - _ Tm 2 f 1 , 



x 



|(l-a) -A{D-l)k 2 a{l-a)-{D-2)ik-u{l-2a)+u 2 + ^ - ^ j 



x exp 



T 

—ik ■ uTa — k 2 Ta(l — a) u 2 

4 



(69) 



Some comments are in order with respect to what has taken place in producing the above 
result. The most obvious one is that a assumes the role of a Feynman parameter as per the 
standard diagrammatic manipulation. Given its introduction the following practical rules of 
substitution, with regard to the original worldline expression, emerge 

G B (r 1 ,r 2 ) = l-2a, G 2 B {r u r 2 ) - G 2 F (n, r 2 ) = -4a(l - a). (70) 

The last relation is especially notable as it explicitly underlines the absence of divergencies 
in the integral over the Feynman parameter. 

Turning our attention to the 7-matrix part, J 7 , of the second order propagator we need 
the following results pertaining to expectation values 

< k ■ ±(T)7 • x(n) >= ^ )D/2 h ■ kG B (r i: T) + k ■ i c/ (T) 7 • x d (r t )}, i = 1,2, (72) 
< x(T) ■ x(n) >= J [JG B ( r< , T) + x d (T) ■ i c '(r,)], i = 1, 2 (73) 

and 

< 1 -x(T)x d (r l ).x d (r 2 > = ^ )D/2 [^7 • x d (T)G B (r 1} r 2 ) + 7 • x^Gb^, T) 

+7 • x d [r 2 )G B {j u T) + 7 • x c1 (T)x c \ti)) ■ x d (r 2 )]. (74) 

As before, we must integrate by parts and make the variable change which introduces 
the Feynman parameter a. Doing so we are led to the following result, in momentum space, 

I-y(q) = ig 2 CF I — wjT^ / dTT 2 f da(l — a)exp[— T{m 2 + q 2 — 2q ■ k + ak 2 )] 

J ( Z7T I K JO JO 

x{ 7 • q[(D + 2)m 2 + (D - 2)q 2 - 2(D — 2)q-k]+j- k(D - 2)(m 2 + g 2 )}.(75) 
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Combining (69) and (75) we finally obtain 



G® 



-g 2 c F - 



d D k 1 



:{m[Dm 2 + (D - A)q 2 - 2(D - 2)q ■ k] 



(q 2 + m 2 ) 2 J (2ir) D k 2 (q - k) 2 + m' 
-if ■ q[(D + 2)m 2 + (D — 2)q 2 - 2(D - 2)q ■ k] - i 7 • k(D - 2)(m 2 + q 2 )} 



m — 17 • q 



-9 2 Cf 



I 



d D k 1 m 



— 2(7 • q — 7 • k) 
{q - k) 2 + m 2 7m 



1 



m — vy ■ q 



(76) 



which explicitly shows that the worldline perturbative, to second order, term for the fermionic 
propagator coincides with the one obtained via the Feynman diagrammatic expansion. 



7. One loop calculation in the fermionic propagator 



In this section we shall apply the worldline approach to a loop calculation pertaining to 
an open line which extends a recent closed line result by Schubert [27], see also Ref [28]. 
Specifically, by going back to (35) we shall perform the loop integration, equivalently, the 
integration over d D k, thereby extracting an expression for the second order contribution to 
the fermionic propagator in parametric form. So as not to bore the reader with endless 
manipulations we shall restrict our attention to the mass term entering (35) which reads 



-9 2 c F 



r d D k 1 r°° m r T r T r 
/ 7^75 12/ dTe' Tm2 drj dr 2 9(r 2 - n ) L )= „ Vx{r) 
J \2%) u Ar Jo Jo Jo Jx(t\=i, 



x{T)=y 



x{x(t 2 ) ■ x(n) - GUn, t 2 ){D - l)k 2 }exp 



1 f T 

-- / drx 2 - ik ■ (x(ti) - x{r 2 )) 
4 Jo 



(.77) 



Introducing a small mass A for infrared protection, as well as a mass [i for the purposes 
of dimensional regularization we perform the /c-integrations making use of the formulae 



A-D 



d°k (k 2 ) 
(27T) D k 2 + X 2 ' 



—ik-x 



A-D 



(4n) D / 2 Jo 



dLL^e" x2L \^ 



D 
2L 



x 

II 2 



e 4L , 



{7i 



where the parentheses around k 2 serve to point out that for numerator equal to 1 the terms 
in parentheses inside the integral on the right hand side are also substituted by 1. 
Isolating the path integral component of I m , i.e. 



1 = L(o)=x Vx(t) 

Jx(T)=ii 



>x(T)=y 

xexp 



x(t 2 ) ■ A(ti) + 



D 



AL 2 



-G 2 f {t u t 2 ){x{t 2 ) -x^ f 



(D-1)D 
2L 



G 2 f (ti,t 2 ) 



1 r T 1 

~4 Jo dT±2 ~ 4L^ Tl ) ~ X<yT2 ^ 



(79) 



we find ourselves having to calculate expectation values of bosonic super-particle field coor- 
dinates with respect to the action 



S = \ f drx 2 + I(x(r 2 ) - x{n))< 



(80) 
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Leaving, for what follows, a number of technical details to the Appendix we present the 
equations of motion in the form 



= da'B(a',a)x*(a'), (81) 

where 

B(a u a 2 ) = [Sfo - n) - 5(*! - r 2 )][5(a 2 - n) - 5(a 2 - r 2 \. (82) 
Next, we introduce the Green's function A^(a± : a 2 ) by 

/ datffa - a)— - -B{a, (tOJAW (<j,<j 2 ) = -Sfa - a 2 ), (83) 
Jo aa A L 

obeying the boundary conditions A^(0,a) = A^'(T,(j) = 0. 

The solutions to the classical equations of motion are readily determined as 

= ~ j(r 2 - r 1 )^[A( 1 )(a,r 2 ) - A^fon)], (84) 

with r/ M (cr) as given in (56). 

The corresponding expression for the classical action is 

Sd = T« • ^ = ^ + 7" 2 7 ( V/" l)2 v ( 85 ) 
4 v 7 4 4 L + G B (7-i,r 2 ) V ' 

Making the substitution x — > x + x cl we have that 

< i(ri)x(r 2 ) >^< x(ri)x(r 2 ) > c +i d (ri)x c/ (r 2 ) 

< (x(r 2 ) - x(r0) 2 >-< (x(t 2 ) - x( n )) 2 > +(x c/ (r 2 ) - x d (r 1 )) 2 : (86) 

where 

< A > = M' 1 [ Vx(r)Ae- s[x] (87) 

J x(0)=(x(T)=0 

with 5 [a;] as in (80) and M~l the usual renormalization factor. 
From (84) we readily determine 

x d (n) = x d (r 2 ) =u, + I« M ( Tl - r 2 ) ^ (ri ' T2) (88) 

and 

[^(r 2 ) - ^(n)] 2 = (r 2 - n) V f — -A ^ . (89) 

In the Appendix we establish that 

AT- 1 < x{r,)x{r 2 ) > = DG B (n, r 2 ) - ^ y^T^ i ( 90 ) 

2 L + G B (ri, t 2 J 
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and 



Af < (x(r 2 ) - x(ti)) 2 > = 2D- 



LG b (ti,t 2 ) 



(91) 



L + G B (r u r 2 y 

where M is defined in the Appendix. 

Putting everything together we obtain (we suppress the Ti,t 2 of the G"s for notational 
economy) 



-mg 



2 °F 



(Air) 



D 



poo n poo , n pT rT 

/i 4 " D / dLe~ x2L / dTT~ D ' 2 e- Tm2 / d n / dr- 

JO JO JO JO 



0(r 2 - n) 



{DG B ° Gl 



2 L + G B 
L 2 



1 + 2 iT2 - Tl) LTG- B 



+ D ~ 1 G 2 



2dJ^4-+ 



L + G b 



+ (r 2 - n)V 



(L + Gi 



2L ° F 



exp 



Tu 2 



u 2 (r 2 - n) 2 
4 L + 



(92) 



After an integration by parts, which rids us of the presence of Gb(ti, r 2 ) in favor of the 
combination — G% we get, once performing the substitution r 2 — Ti = aT, T\ + r 2 = 
2/3T, L -> LT, 



(4^ L 



dL / dTT 3 - D e- T ^ m2+XL ^ 



da- 



1 — a 



L + a(l — a)\ D / 2 



x 



u 2 (D - 1) 



cr(l — a) 



[L + a(l -a)} 2 



2D(D-1) a{l-a) (D - 2)u 2 a{\ - 2a) 



T L + a(l-a) 



D 1 - 2a 
T 1-a 



exp 



Tu 2 Tu 2 



a 



4 4 L + a(l-a 
The corresponding expression for the Fourier transform reads 



2 L + a(l-a) 
(93) 



(4vr)^ Jo ""Jo Jo ^(L + a) D / 2 

x [D(m 2 -q 2 + \ 2 L) - 2(D - 2) Lq 2 ]exp[-T (m 2 + \ 2 L - q 2 



a 



L + a 



= —mg 



2 " (1 da r D S m2 7 a2) ^ D : m -«}ts{2-D/2 )m 



(47T 



/i 



o [m 2 a + (1 - a)X 2 + q 2 a(l - a)} 2 ~ D l 2 



As already mentioned, we shall refrain from working out the gamma-matrix part of the 
propagator. The point is that our methods have produced a one loop result in a parametrized 
form pertaining to a configuration with an open fermionic line, thereby extending the scope 
and use of the approach and methodology of Ref . [7] . 



8. Concluding remarks 

Our efforts, in this work, have centered around the role of the spin factor in the worldline 
approach to gauge theories with fermionic matter fields. Placing specific emphasis on open 
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particle path configurations we were able to accommodate the spin attribute of the prop- 
agating entity into a geometrical type of expression (spin factor) which registers as a well 
defined weight factor in the path integral. Further to the point, its presence facilitates the 
extraction of simple computational rules which, at the perturbative level at least, produce 
algorithms of the Bern-Kosower type [14] while extending, at the same time, their applica- 
bility to fermionic Green's functions. On the other hand, the dynamics of the particle-based 
casting of the system amounts to emission and absorption of gauge field modes by the matter 
particle which is solely controlled by a Wilson line operator. Our scheme thereby achieves a 
neat separation between geometrical and dynamical aspects of physical descriptions in gauge 
field theories. To summarize, our formulation has not only made it possible to accommodate 
open fermionic paths, equivalently, Green's functions, into worldline considerations but has 
produced computational rules that have clear as well as simple interpretations. 

The fact that our approach has identified the Wilson line operator as the sole agent of 
the dynamics operating in the system is worth some further comments. This occurrence has 
proven to be of central importance to the dynamical considerations we have been pursuing 
in connection to long distance properties of QCD [10-12]. The wider realization that Wilson 
line operators carry the dynamics of eikonally-based considerations in QCD is, actually, well 
founded in a variety of approaches which aim at factorizing soft contributions to quantities 
of physical interest, such as form factors and parton-parton scattering amplitudes [29-34]. A 
marked contrast is that in the worldline approach to QCD Wilson lines enter directly as an 
indigenous component of the action whereas in standard field theoretical approaches they 
are brought in, within the context of an operator formalism, for the purpose of representing 
quarks as ordered exponentials of the gluon field connected at the point where the hard 
scattering occurs. In other words, they enter as sort of 'tails' attached to the quark field 
operator for the purpose of accounting for long distance behavior. Our present work further 
establishes the role of the Wilson line operators as the quantity which carries the dynamics 
in a gauge field theory even beyond soft effects. 

Be that as it may, our immediate plans are to study situations wherein external gluonic 
fields attach themselves to an open fermionic line thereby extending the present considera- 
tions to vertex functions and 'Compton' amplitudes. We expect to report on progress along 
these lines soon. 
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Appendix 



We display, below, some interim steps associated with the derivation of formulas given in 
Section 7. Consider the particle-based action functional S as registered in eq. (79). Casting 
it in the form 

1 f T 1 f T f T 

S = - dTx 2 (T) + — do x \ da 2 B(a 1 ,a 2 \r 1 ,T 2 )x(a 1 ) ■ x(a 2 ), (A.l) 
4 Jo 4L Jo Jo 

with B((Ji, <t 2 |ti, t 2 ) as given by eq. (82), we are led to the equations of motion given by eq. 
(81) in the text. 

Considering a solution of the form x^ — rj^ + z^, with rj^ as in (56) we are led to the 
relation (we drop, as in the text, the r-arguments in B) 



/ da' 
Jo 



u. 



z,(a') = [5(a - r 2 ) - S(a - n)}^ - r 2 )-£ (A.2) 



which justifies the introduction of A^\a, a') according to (83). One, now, readily determines 
=Vv~ [ da'A^{a,a')[5{a' - r 2 ) - 5{a' - n)](n - r 2 )^ (A.3) 

which produces eq. (84) in the text. 

Employing the notation 5(a — a')-^- — jB(a' , a) =< a\(d 2 — j^)\a' > we find, following 
the procedure applied in Refs [27,28] for the closed contour case, 

-l 



-A^(a,a') =< a\id 2 - 



B 



W >=< a\ 1 - 



B 



d~ 2 W> 



=< a\d- 2 W >+-< a\d- 2 Bd- 2 W > +— < a\d' 2 Bd' 2 Bd- 2 \a' > +... (A.4) 

with the first term in the series to be denoted, for the present purposes, by A^(a,a'); it 
actually coincides with the free bosonic propagator, cf. eq. (50). 

Straight forward manipulations on the matrix elements entering the above series result 
to the following, closed-form, relation 

A<-W) = A'°W) - |A1 °' ( "- ) - AW '^»f'^ - (A.5) 

L + Gb{ti, t 2 ) 

Liberal use of (A.3) and (A.5) leads to the relation 



L + G b (t u t 2 ) 



(A.6) 



which, in turn, leads to eq. (89) in the text. 

From (A.3) we obtain, after some straight forward manipulations, for the derivative of 
the classical bosonic particle fields 



^tW) =u ll - -iv(n-T 2 ) 



G^jayri) - G B (o r ,r 2 ) 
L + G b (t 1: t 2 ) 



(A.7) 
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which justifies eq. (88). 

Turning our attention to expectation values, let us first furnish an outline derivation of 
(90). Under the substitution ni: d + iwe have that 



< x(n) ■ x(r 2 ) >= [< i(ri) • x(r 2 ) >o +x c1 (t 1 ) ■ x c \t 2 ) < 1 > }e 5C \ 



where 



< 1 > 



Jx 



x(0)=x(T)=0 



Vx(r)exp 



£ drx* - ±(x(t 2 ) - x(n)Y 



l r T 

4 



One easily obtains 



< 1 > 



° (4vrT) D / 2 



L + G B (n,T 2 ) 



D/2 



(A.8) 
(A.9) 

(A.10) 



On the other hand 

< x(ti) ■ x(t 2 ) > = 

xexp 



Vx(t) 



SJ^(ti)5J^(t 2 ) J x (o)=x(t)=o 

I rT I rT rT rT 

- \ drx 2 - — / do x \ da 2 B(ai, a 2 ) - i / drJ( 
4 Jo 4L Jo jo jo 



T • X(T) 



(A.ll) 

J=0 



Now, from (A. 5) we determine 

d 2 
dada 



A {) M = - 2 GbM-- a L + G B (n,T 2 ) • (A ' 12) 



which facilitates the carrying out of the (functional) Gaussian integration in (A.ll). One 
finds 



< x(ti) ■ x(t 2 ) > 



1 



D/2 



(4ttT) b / 2 {L + Gsin,^ 
xD 



fi ( _ \ 1 G l( r l' r 2) 

(*b{Ti, t 2 ) - -- 



2L + G b (t 1 ,t 2 ) 



(A.13) 



as per eq. (90) in the text. 

Finally, eq. (91) is deduced as follows 

< (x(t 2 ) - x(n) 2 > Q = < (x 2 (r 2 ) > + < (x 2 (n) > Q -2 < (x(t 2 ) ■ (x(n) > Q 

= 2D[A« (r 2 ,r 2 ) + ^(n.n) - 2A«(r 1 , r 2 )] < 1 > 

- 2D L + G b (t 1 ,t 2 ) <1> °- 



(A. 14) 
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